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On Interpolation ; an Essay containing a simple Exposition of the 
Theory in its most useful practical applications, together with a 
general and complete Demonstration of the Methods, of Quinti- 
section of Briggs and Mouton for equal intervals ; and of the 
process explained by Newton in his " Principia," for intervals 
of any magnitude whatever.* By M. Fred. Maurice. 

Translated from the French by Thomas B. Sprague, M.A., 
and J. Hill Williams, Vice-Presidents of the Institute of 
Actuaries. 

1 HE calculation of extensive tables of logarithms could not fail to 
give rise to the method of interpolation ; and such indeed was its 
origin. After Napier had made his memorable discovery, and 
he and Briggs had recognized the advantages that would result 
from the adoption of 10 as the base of a system of logarithms, it 
was Briggs who courageously undertook the immense labour required 
for calculating those famous tables, published in London in 1624, 
which have not been surpassed in extent and accuracy by any sub- 
sequent work. It is well known, however, that Briggs left a great 
hiatus in those tables, and that the logarithms of numbers between 
20,000 and 90,000 are not to be found in them ; but scarcely were 
the first tables printed, when Briggs undertook with fresh energy 
the computation of the Logarithms of the Trigonometrical Lines, 

* From the Connaissancc des Temps for the year 1847. 
VOL. XIV. B 



2 On Interpolation. [Oct. 

and he was on the point of terminating this vast enterprise when 
arrested by death. It was his friend Gellibrand who completed 
the work, and gave it to the world in 1633. 

The great service rendered to science by this distinguished 
computer is not confined to the exact determination to a great 
number of decimal places of the logarithms calculated by him, and 
of their differences of several orders. Briggs has besides shown 
himself to be a mathematician of great skill, by his learned pre- 
faces to the two great works to which we have referred. In these 
he has developed a host of ingenious methods invented by him- 
self, with the view of combining rigorous exactness with facility of 
computation ; and, with that object, he has given various remark- 
able methods which enabled him to verify his results to any degree 
of accuracy required. 

One of the most remarkable of these is the singular process he 
explains(chapter xiii. of the Preface to his first tables, and chapter xii. 
of the Preface to his Trigonometrical Tables) for interpolating four 
intermediate values between each adjacent two of a series obtained 
by direct calculation, both the original and the interpolated values 
being separated by equal intervals. To this explanation he adds 
that, in publishing to 14 decimal places the logarithms of the 
30,000 numbers comprised between 1 and 20,000, and between 
90,000 and 100,000, it was his intention to give, by this process, 
an easy means of determining the 70,000 logarithms remaining to 
be calculated in order to complete those of the first hundred 
thousand natural numbers; and it may well be thought that he 
had already employed it to compute some of the logarithms con- 
tained in his valuable tables — we mean the Tables of Logarithms 
of Numbers, edited by himself in 1624. 

This method, which is very remarkable, and preceded by nearly 
half a century the researches of Mouton and Regnauld upon a 
similar method of interpolation, is given by the author without 
demonstration. He only lays down the highly complicated rules 
of his method, and that in a manner wholly wanting in the clear- 
ness which, since the time of Euler, analysis has learnt how to 
throw on most subjects of which it treats. Furthermore, the 
process itself is more singular than easy of application : there 
seems nothing in it to guide the computer in choosing the limit of 
the corrections to which he must confine himself; and, moreover, 
the differences, of various orders, of the function under consideration, 
having reference to different values, might give rise to some con- 
fusion in practice. 
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To this defect in the explanation we must, no doubt, attribute 
the silence of mathematicians as to an invention so remarkable and 
of such merit. No one drew attention to its originality until, after 
the lapse of two centuries, Legendre made it the subject of a learned 
paper in the additions to the Connaissance des Temps for 1817. 
He there demonstrated, by means of a rigorous modern analysis, 
the reasons of the rules given by Briggs for deducing the differ- 
ences to be employed in interpolation from those which he calls 
mean differences ; showing that these last have, in fact, necessary 
relations with the differences of the original terms. 

In that paper, the shortness of which is a matter of regret, 
Legendre confesses that his own demonstration is wanting in sim- 
plicity, and expresses a hope that some other proof may be found 
more closely akin to that which the author must undoubtedly have 
discovered, although he did not publish it. In the present paper 
we purpose to answer that appeal ; and we believe we have succeeded, 
without using any process or method of investigation not known to 
mathematicians living in 1620, contemporary with Kepler and Har- 
riot, and before Descartes and Fermat. 

With this view, it appeared advisable, before giving our 
elementary demonstration of Briggs's curious process, to examine 
all the methods of investigation which were within his reach when 
about half way through his great work ; and this will be the object 
of our first section. In the second, we shall describe, in Briggs's 
own words, that process which Legendre has deemed so worthy of 
attention ; and in the third, it will be completely demonstrated, by 
a purely arithmetical method, whose only resemblance to modern 
analysis is due to the convenient notation so generally used at the 
present time. 

We have contrasted with this rather obscure method the clear 
and convenient one published in 1670 by Mouton. He was a 
priest of Lyons, and the author of an astronomical work on the 
diameters of the sun and moon, containing several ideas which, con- 
sidering the period in which he wrote, are very remarkable. But 
this method, which Mouton states was generalized by his friend and 
countryman Regnauld, is only described by him at some length, 
and is not demonstrated. Lalande was the first who undertook (in 
1761) to demonstrate it for the three first orders of differences; 
and no further progress had been made, when Lagrange, in the 
Berlin Memoirs for 1792-3, published a demonstration as general 
as it was learned. Shortly after, Prony took this excellent method 
for the basis of the immense work executed under his direction, for 
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computing the great logarithmic tables called Tables Du Cadastre 
(Government survey), and, we are told, developed it in the preface in 
all its details. Those tables, however, which surpass both in extent 
and accuracy all that have been published, only exist in two manu- 
script copies, carefully preserved in different places, so that l'rony's 
analysis has never been given to the public. It is true that Lacroix, 
in the third volume of his great treatise, has briefly given an outline 
of the method and of its demonstration ; but as it is far from elemen- 
tary, we trust that our readers will not consider as useless the more 
simple and complete exposition of it which we give in the fourth 
section. We shall there further show that Briggs's and Mouton's 
methods, notwithstanding their great dissimilarity, lead to exactly 
the same results, not only in their numerical applications, but like- 
wise in the general formulae. As, moreover, we have not omitted 
to give such details as may be practically useful in interpolating, it 
will be seen that the clear and easy process of Mouton offers unusual 
advantages from the symmetry of all its operations. 

Hitherto we have only considered the case when the intervals 
between the given values are equal ; but other cases occur in prac- 
tice. For example, in some astronomical observations it happens 
that the times between the successive observations are unequal ; and 
then, taking the time as the variable, the values given by observation 
are not equidistant, and the problem is to find from the given data 
the value of the observed element at any given time. 

Newton gave the first solution of this problem in a lemma in 
his Principia, which, according to his custom, he did not stop to 
demonstrate; and Laplace, who adopted that solution in his method 
of calculating the orbit of a comet, also confined himself to its 
simple enunciation. We devote our fifth and last section to the 
general consideration of cases of this kind; and after giving a 
solution based on the same purely algebraical reasoning we have 
used throughout, we shall conclude with an easy demonstration of 
Newton's solution, under the form given to it by Laplace in his 
Mecanique Celeste. 

Lastly, in an Appendix, we have considered more particuhrrly 
the formulae most frequently used in astronomical interpolations; 
and we have deduced them also from the same simple principles. 
We would particularly draw attention to the manner in which we 
there treat a formula which Stirling, in his excellent work on 
interpolation, had only obtained by induction, and to the additions 
we have given to M. Bessel's demonstration of that useful formula. 
The combined simplicity and clearness of all these calculations 
will, it is hoped, secure for them a careful perusal. 
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§1. 

1. Some persons who have carefully examined the prefaces of 
Briggs, have seen therein reason for thinking that he had antici- 
pated Pascal, Wallis and Newton, in the discovery of the celebrated 
Binomial Theorem. Dr. Charles Hutton, for instance, in the 
introduction to his Tables, states this explicitly, and endeavours 
to explain how the fact may be reconciled with Newton's acknow- 
ledged claim to that most valuable discovery. We shall see, how- 
ever, that Briggs may have arrived at a very analogous formula by 
a widely different course of reasoning, which would not necessarily 
require him to cousider the raising of a binomial to a given power, 
either integral or fractional. We shall do our author full justice 
by allowing that he was led by the very nature of his investiga- 
tions to find general formulae showing the relations between the 
equidistant terms of any given series and their differences of any 
order. It is evident that the logarithms of successive numbers 
are quantities of this nature. 

2. Thus, if we have a scries of given quantities 

ti u , «,, « 2 , it 3 , ....,«,,..., 

corresponding to values of the variable, x, 

0, 1, J, o x, . . . , 

it is easily seen that we may assume the general term of the series 
to be given by the formula 

u s =:ti + Ax+B.x(x-l) + C.x(x-l)(x-2) + T).r(x-l)(x-2)(x-3)+ . . . 

where the coefficients A, B, C, D, . . . . are to be determined by 
means of the conditions (obtained by putting x equal to 1, 2, 
3, . . . . successively) 

« 2 =w + 2A + 2B, ? < 4 =» + 4A+]2B + 24C + 24D, , 

Hence we deduce, representing by the symbol 8 the difference of 
two of the quantities, and by 8 n , generally, a difference of the nth 
order, 

A:= «[— v = Z>i ; 

« 2 — 2Ui+U V 2 — ?/, — (J/, — tl a ) <S(W| — u n ) 2V 



B= 

c 
I) 



2 2 2 1.2 

"3— 3"2 + 3"i— "n t 3 "n 



1.2.3 1.2.3' 

w 4 — 4»3+C;/ 2 — 4« l +;?| l _ < A it n 
1 2.3.4 ~ 1.2.3.4 
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Substituting these values, the formula for the general term becomes 

,,,x x * x(x—l) x(x—l)(x — 2) „ 

x(x-l)(x-2)( x-3) 
+ iJ^A ^o+..-. 

3. Tins expression calls for two remarks : — 

The first is, that the series may coine to an end, either because ,r is 

a finite intecer.and then all the coefficients after — .—r— ^ 

1. .2 ... x 

will of necessity vanish ; or because we meet with constant 

differences of the ?jth order, in which case all the terms after 

— . ' ' ' • S"(/ will disappear, because g n+, M = O, &c. 

The second remark is, that we have not proved that the form 
which we have seen to prevail in the case of a few terms, will be 
generally true whatever the number of terms. But a very simple 
consideration will satisfy us on this point. Suppose that x is 
increased by unity ; we shall have 

x . x(x— 1) ,, x(x— \)(x — 2) „ 
■•• «x+i=»o+ j • K+ J Y2~ * + i 23 ' 0+ • • ■ 

. i * x *•> x(x—l) 

+ l.cw,, + - • e-u + l 2 • t 3 "o + • ■ • 

x+l „ (x+l)x ,„ (r-i-l)r(j-l) „ 
= '<n+ — —■"<<)+ —J2 C "" oH TJ3 •c 3 u + ... 

Thus u x+1 is expressed by a scries in perfect conformity with 
that of the formula (E) ; from which it may be strictly concluded 
that if the formula is proved for any term, u x , it is true also for 
the next term, v<. +1 ; and by successive inductions it is true for all 
positive and integral values of u: 

This formula (E) is so related to that of the binomial theorem, 
that a century and a half after Briggs's time it was at last written 
in the abridged form 

U £ ={l+h}'l t0 . 

In fact, the formula for v x results rigorously from the product of 
w by the expansion of the ,rth power of the binomial (1 + §). 

4. We are able to see also, when the coefficients A, B, C, D, . . . 
have been determined, that the successive differences of u are 
expressed in series containing the successive values of u itself, 
with alternate signs and coefficients similar to those of the series 
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(E) ; so that, for example, up to the nth order inclusive, we have 

obtained the expression 

/T% . . n n(n — l) w(«— 1)(» — 2) 

(D) S"u = «„- - ■ «„_, + ^j-^-J • «„_ 2 L2 A 3 -' • «_, + . . . 

A line of reasoning similar to that of the preceding article will 
prove that this general form is also true when n becomes n + 1. 

In fact, S B + , « = S B (« 1 — «o). But the value of S"w, will be 
deduced from that of S"« just found, by writing u n+l for u„, w„for 
«„-!, &c. We shall then have 



n n(n— 1) «(»— 1)(« — 2) 



d" +1 «/ =!«„ +1 — - • u„+ — p^— • «„_i 7-3-3 «.-»+ • 



n n{n— 1) 

ra+1 (n+.l)» (fi + l)«(»-l) 

But this expression being of exactly the same form as (D), we 
see that if that formula is true for the rath difference, it is also 
true for the (w+l)th difference; and by successive inductions it 
is true for any order of differences whatever. 

This new formula received, about the same time as the previous 
one, the following concise and convenient form — 

2"m=(m-1)", 
where care must be taken, in the expansion of the second member, 
to write throughout u n instead of u K , u n _ x instead of w" -1 , &c; 
and in the last term, instead of 1 or u°, to write u . 

5. We have considered, (art. 2) a series of given quantities, 
« , M,, « 2 , . . . , u x , . . . or values of the function u x , corresponding 
to the values of the variable 0, 1, 2, . . . , x, . . . ; and we have 
obtained the formula (E) for the general term u x . 

Let us now consider the series u , u,,, uy t , ...,«*,... where h 
is an integer, and where we must have z=xh, x being also an integer. 

Then it is evident that the formula for the general term u, will 
be deduced from that for u x or from (E) by substituting for x in 

the second member its value j; for x takes the values 0, 1, 2, . . . 

x, . . . when z is according to our hypothesis made equal to 
0, h, 2h, . . . xh, . . . successively ; only it will be necessary to put 
Aw or u h —u , and its successive differences A 2 tt > A 3 m , ... in the 
place of Sw , or «, — u 0> and its successive differences &u , 8 3 m , . . . 
Making these changes, the formula becomes 
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._. z z(z — A) .„ z(z — A)(z — 2A) , 

which may also be obtained directly, in just the same manner as 
the formula (E) was first obtained (art. 2), if we represent by Aw 
the difference u h — u . 

6. This formula (F), which for values of z equal to 0, A, 2A , . . 
gives 

"* =«o + A« 

%*=%+2A%+A 2 « 

&c. = &c. 
will of necessity coincide with that deduced from the formula (E), 
by making x equal to 0, A, 2A, . . . successively. Conversely, 
(E) must also coincide with (F), when z is made equal to 0, 1, 2, . . . 
successively, which gives 

t "o="o 

1 A 1(1 -A) A0 1(1-A)(1-2A) , 



(zy J 



2 A 2(2— A) , 2(2-A)(2-2A) ; 

. W + A A " 0+ T2JT A2 " 0+ 1.2.3.A* A3 »° + • • • • 
&c.= &c. 

We thus conclude that tbe use of the single formula (F) will 
give the values of v 2 for all integral values of z. Consequently, if 
we have computed a series of values of a function u z , such as are 
shown in (Z), we see from (Z)' that we can determine, by means of 
the known quantities, u„, Aw , A 2 w .... A x_1 w , the A — 1 values 
of u z comprised between u and u h , or between u h and « 2 a. &c-> and 
this is called interpolating ; i.e., if we have a series of values of u z 
corresponding to values of z, that differ by the constant quantity A, 
we know how to insert, between any two adjacent values, (A — 1) 
intermediate values, equidistant like the first. 

7. Similarly if the original successive values of z had been 

0, 1, 2, 3, ... , instead of 0, A, 2A, ... we should have had a 

formula exactly similar to (E), and which would be written thus : — 

/,^x z* s(z— IKo z(z— l)fs— 2)„ 

(E') u,=u n + T 3u + -i_2 ^ Mo + -^— ^ i 3»i/ + • • • • 

and the known values of u, would be 
«o = "0 

«1 = «0+$"o 

"2 = «o + 2£«o + c 2 Wo 
&C.= &C. 
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But if we now suppose z to have the values 0, y, -j, . . . -= , . . . 
successively, it will be evidently sufficient to put, for z, in (E') its 

X 

general value yj for (E') is really the same formula as (E) : they 

are both dependent upon increments of the independent variable 
equal to unity ; so that in the form which (E') then takes, incre- 
ments of x equal to unity will give the successive values of u z . 
We shall have then 

,-n,-. x. x(x — h) .„ x(x — h)(x — 2h) „„ 

( F<) w =„ 0+ - *,+ ±~j ax+ -^y- **+ • • • • 

a formula, which, making x equal to 0, h, 2h, . . . , will give again 
the preceding values of « , «,, w 2 , ■ . . , namely, 

Consequently, if we have a table showing, like the formula (E) 
the values of u x for values of x equal to 0, 1, 2, 3, ... , so that we 
know the values of u , Su 0) S 2 u , S 3 u , . . . , the formula (F') will 
give the {It — 1) values mi, W2, . . . tw-i, intermediate between w 

h h h 

and Uh or u x , or enable us to interpolate between u and u x these 

[h — 1) values which are, like the original values, equidistant. It 
is evident that we may do the same for the second interval com- 
prised between w* or Wj , and un or u 2 ; for the third, and for all 

ft h 

subsequent intervals. 

8. Such are the means of calculation which we may suppose 
Briggs to have possessed when he began to work at his great 
tables. In these operations we find really nothing beyond the 
general arithmetic then known ; as far as theory goes, they only 
introduce the simple and precise ideas of the differences of various 
orders of given quantities, and the reciprocal relations between such 
quantities and differences. 

We shall, however, soon see that these formulae and ideas were 
more than sufficient to lead him to the singular method of inter- 
polation which we are about to explain, and which we shall first 
present as it appears in the 13th chapter of the Preface to the 
Arithmetica Logarithmica, published at London in 1624. 

9. We have said that these means were more than sufficient to 
lead Briggs to his method. We shall see in fact that it is highly 
probable that he never obtained the formulas (F) and (F') ; at ail 
events, they are not required in order to obtain his results. More- 
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over, it may be remarked, that if he had only had before him the 
second of the formulae (Z') in art. 6, which gives 

1 1(1 -A) „ 

he would have been led to think that S"(/ must depend upon A"u , 
A" +1 (/ > •••> an( l that, consequently, if A"u were a constant 
difference, it was unnecessary to notice S" +1 « > S" +2 tfo> • • • • But, 
as will be seen in articles 11, 18 and 28, Briggs carries the formula 
for <$'</ up to the ten thousandth part of S 20 w , . . . , whilst, sup- 
posing A 5 M constant, he should only have taken S r '« into con- 
sideration, and should have neglected all differences of a higher 
order. 

But these observations will arise more naturally as we proceed. 
Let us now hear Briggs's own description of his curious process. 

§11. 

10. In the 13th chapter of the Preface to his Tables, Briggs 
wished to give the means of finding the logarithms which he had 
not calculated, namely, those of the numbers between 20,000 and 
90,000 ; and for that purpose he proposed the following problem : — 
" Given a scries of equidistant numbers and their logarithms, to 
find the logarithms of the four numbers interpolated at equal 
intervals between each adjacent two of the given numbers." 

Thus, for example, Briggs having calculated the logarithms of 
all the numbers below 20,000, easily obtained, by adding the 
logarithm of 5 to those of the numbers 4,001, 4,002, 4,003, . . . , 
the logarithms of the numbers 20,005, 20,010, 20,015, . . . ; but 
he had still to find a means of obtaining, without too much labour, 
the logarithms of the four whole numbers comprised between each 
two of the last numbers. 

But let the author speak for himself. " Take," he goes on 
to say, " the first, second, third, fourth, &c, differences of the loga- 
rithms of 2,115, 2,120, 2,125, . . . , and divide the first differences 
by 5, the second by 5 2 , the third by 5 3 , . . . (or, which comes to 
the same thing, multiply them respectively by -2, by -04, by -008, 
&c.), and these quotients (or products) will be what I call the 
mean differences of the corresponding orders." Briggs then gives 
examples of the manner of obtaining them ; and afterwards pro- 
ceeds to correct them, so as to be able to employ them in his 
investigations. 
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11. Having found that the fifth differences of the logarithms 
of the numbers 2,115, 2,120, 2,125, . . . , are constant, Briggs 
observes that the mean differences of the fourth and fifth orders 
cannot be corrected; for, he says, the sixth and seventh differences are 
nil; adding, " omnis autem differ entiarum correctio fit per subduc- 
tionem differ entiarum alter narum magis remutarum et correcta7-um" 
(every correction of the differences is made by subtracting the alter- 
nate corrected differences of higher orders) (Preface, p. 27, near the 
bottom). This passage appears at first sight to involve a vicious 
circle, inasmuch as it requires, for the calculation of the corrections, 
that we should have differences already corrected. But this diffi- 
culty disappears when we notice that Briggs begins by obtaining 
differences which, requiring no correction, may be said to be 
already corrected ; so that when he subsequently makes all these 
corrections, commencing with the highest order, he obtains suc- 
cessively the corrected differences he requires, whose order in effect 
alternates by increasing at each step by two units (See Articles 15 
to 18 and 28 below, and compare them). 

Briggs contiuues : — " Thus the subtraction of the seventh differ- 
ences corrects the fifth differences, that of the sixth differences 
corrects the fourth, and so on." Therefore, in the case under con- 
sideration, we may take the mean fourth and fifth differences for 
the corrected fourth and fifth differences. 

As to the mean third differences, we are to correct them by the 
subtraction of three corrected fifth differences, and Briggs gives a 
numerical example. 

The mean second differences are corrected by the subtraction of 

7 
two corrected fourth differences, and also of -=■ of the sixth differ- 

5 

ence, if it have any sensible value. Briggs here gives a fresh 

example. 

Lastly, we arc to correct the mean first differences by deducting 

therefrom one corrected third difference, and ■= of the fifth difference: 

5 

and Briggs gives another calculation as an example. 

" Such," adds Briggs, " is the course to be pursued in correcting 

the differences, whatever may be their order, commencing always 

with the highest order." He thereupon gives a table showing all 

the subtractions required to correct the successive mean differences 

of every order from the first to the twentieth. It shows, for 

example, that the difference of the fourth order is corrected by the 

subtraction of the following multiples of the higher corrected dif- 
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ferences, each difference being represented by the exponent of its 
oider placed within brackets : 

r 4(6) + G-8(8) + G-4(10) + 3-64(12)+ 1-28(14) "I 
L + -272(1 o) + -032(18) +-00 10(20) J- 

12. Briggs then gives some practical suggestions, intended to 
prevent confusion in the application of his rules ; and, as an example, 
he finds the logarithms of the eight numbers comprised between 
2,115 and 2,120 and between 2,120 and 2,125. For this purpose 
he forms a table, in which are arranged in different places both the 
given logarithms and their differences of the first four orders, 
corrected as previously explained; those of the fifth order being 
neglected as insignificant. He finally states that the required 
logarithms will be found by adding the fourth differences to the 
third, the sum of these to the second differences, the new sum to 
the first differences, and the sum thus found to the preceding 
logarithm. 

But this table appears complicated and badly arranged, if we may 
venture to say so. For we must subtract, and not add the diffe- 
rences, since the differences of the logarithms of increasing numbers 
decrease in value; and Briggs himself does this. It is, moreover, 
impossible not to observe the needless fulness of a process by which, 
the mean differences of the fourth order being reduced to two 
figures in the thirteenth and fourteenth place of decimals, we are 
told to correct them by means of differences smaller in themselves, 
which, however, have been previously computed to the sixth, eighth, 
.... twentieth, place of decimals. 

It is not our object here to praise this process, which, as 
we have just seen, Briggs greatly modifies in his actual com- 
putations, but only to explain and prove it. For this is all which 
the author has given by way of explanation, and that is but 
little. He only adds that this method, which may be called 
quintisection, applies, with the requisite alterations, to the steps 
necessary for effecting trisection and septisection, but that he greatly 
prefers the method of quintisection. 

§111. 

13. If Briggs had been in possession of the formula (F), he 
might have thought of using it to accomplish the object he had in 
view in his method of quintisection. In fact, making /j = 5, and 

and treating A 5 « as a constant quantity, the quotients— ~ , in that 
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formula would be his mean differences. For brevity, write A" for 

-tj-^j then applying (F) to find the values of u i} u 2 , « 3 , . . . , as 

far as terms involving A 5 , \vc might find Su Q , Su lt <5m 2 > • • • > 8 2 a , 
S 2 «i , . . . , in terms of A 1 , A 2 , . . . , A 5 , whence, conversely, we 
could deduce the values of A 1 , A 2 , . . . , A r> , in terms of <Sw„, S 2 « > • • • 
But there is no reason to believe that Briggs was acquainted with 
that formula; and moreover, in following the course just pointed 
out, he would have encountered long and complicated calculations 
wholly devoid of symmetry. 

14. The formula (D) of art. 4, which is much more simple, 
was more likely to have suggested itself to him, as offering great 
advantages for the object he had in view. In fact, when the 
variable receives successive increments of 5, that formula gives 

[Att = ;/ 6 — ii 

A 2 «o=«io— «5 — («5 — «o) =1> > suppose, 
(A) J jA 3 » ="i5— "io— («io — "s)-P=Q — P 

AX=«2o-"i.'-(«i5-«io)-2Q+P=R-Q-(Q-P) 

A 5 «o=»2o-"20 — («2o— "is)— 3B + 3Q— P 

= S-R-2(R-Q) + (Q-P) 

On the otlier hand, the formula (E) gives, putting ,i' = 5, 
« 5 — tt =5{h( + •2c i it + 2Pti (l + S 4 u + - Pk ), 
which may be reduced to 

(4) » 5 — » =5(hc 2 + Piii + - Pu ). 

But we have, by the first principles of the method of differences, 

w 5 — "0 =&'0 +<>>'l +2» 2 +t>»3 +S«4 
"10 — w 5 =^»5 +«"0 +«'7 +^«8 +2« 9 
«15 — M 10 = ^"l0 + 2«u + 3«] 2 + £«13 + <> U U 
&C. = &C. 

Here it is at once seen that each of the preceding equations 
may be deduced from the preceding one by increasing the variable 
by 5; and we must therefore get expressions for u t0 — m 5 , 
« 15 — u l0 , . . . precisely similar to the value given by (b) for 
w 5 — w . Thus, 

M io— «5 =o(Sn 7 +Pu 6 + -S 5 m 5 ) 
o 

«i& - u l0 =5(Su l2 + 2 3 «n + - S s «io) 
&c. — &c. 
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But we can readily satisfy ourselves that we shall have likewise 

"7 — !' 2 = 5((5" 4 +<5 3 « 3 + 7<J : '!/ 2 ) 
«'13— «8= 5 ( 2 "l0 + ^ 3 «9+ 7 C 5 (/ 8 ) 

&c. = &c. 
and we shall be justified in concluding generally that 

(i5) «».— «m-5=5(2»,„_ 3 +a 3 i/„,_ 4 + -a 3 « m _ 5 ) 

from which we can easily deduce 

(C) ^(«m-"».-6) = 5(3" +1 "„-3+2" +3 «»-4+ Jc"+^ m _ 5 ). 



15. These preliminaries being settled, let us return to the 
relation (b). It will give, dividing by 5, 

— = — = -=- =A 1 =dM 2 + d 3 « 1 + - e-'(/ ; 
5 5 y 

from which we obtain 

hi 2= A 1 — C 3 «! — - S 5 M - 

o 
But this is the first expression in Briggs's table, which has been 
verified by Legendre. 

(16.) The second of the formula? (A) and the relation (B) give 
in turn 

0) A 2 « =5 ■d(u 7 -u i ) + P(u 6 -u l )+ -&(u s -u )\. 

But by means of the relation (C), we get 

1 

5 



2(« 7 -« 2 )=5(a 2 « 4 +S' 1 «3+ i-2 6 « 2 



3 (« 6 -M 4 )=5 ( a 4 «3 + 6 H 2 + - &lh 

l s»(«» - « )=5 ^ J a 6 « 2 + i «"«, + Ji S'X 

Adding and multiplying by 5, 

A 2 « =5 2 (c 2 m 4 + 2S 4 k 3 + 1 -4a c « 2 + -42% + -04a 10 M ) ; 

or, dividing by 5 2 and transposing, 

(2) 3 2 (/ 4 =A 2 -22 4 « 3 — l-4a«« 2 — •42 3 « 1 --04S 1 %. 

This is the second formula in Briggs's table, which has also been 
verified by Legendre. 
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17. The third of the formulae (A) and the relation (B) will 
give similarly 

A 3 »„=5 |S(W|2— «7>+S 3 ( w ll-«6)+ jjS 5 («10 — « 6 )|— A 2 "o! 

or, by virtue of (C) and substituting for A 2 w its value as given 

by(.r), 

(y) A»«o = 5» j««(«'9-«4) + 2 4 («8-«8)+ ^ fi («7-«2) j 

V t) t) «> J 

=52{a2(« 9 - W4 ) + 2av8-"3)+^ («7-«2) 

+ |g 8 ( M6 _ ttl)+ l_a'V 5 -«o) 

=5 3 [s'» + 3 5 m § + J c 7 (/ 4 + 2^ s «, + Pu 4 + J a»«, 



+ ya ,1 «,+« 1, «i+Ja i5 «o)] 

whence, dividing by 5 2 , reducing and transposing, we get 

(3) Pu 6 =A 3 -S¥u 5 -3-6d ,7 u i -2-2Pu 3 --72Pht 2 --12c l3 u l --00Sd l hi , 

and we thus have the third formula in Briggs's table, the last which 
has been verified and demonstrated by Legendre. 

18. We pass now to the fourth of the formulae (A). It gives 

A 4 w = { "20— «is — («n — «io) } — A 3 w — ( A 3 «/ — A 2 «< ) ; 

whence, by means of (B) and (C), substituting the values of A 2 w 
and A 3 w , as found from (a?) and (»/), 

A 4 K = 5 | S(«|J — «ia) + t> 3 (!<iG- «i l) + T 5 0l5 - «1<>) 

— 3(m 12 — « 7 ) — 2 3 ( «„—«,.)— -#>(«,„ — M 6 )L 

— 5|(X«12 — «7)+0 3 («ll — "«)+ g 2 5 («, — «,) 

— S(«7 — M 2 ) — 1> 3 ("6 — «l)— 7 ^(W,— «o)| 
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= 2 |^(«14-«9) + f 4 ("l3-"8)+5^'("l2 — «7) 

-C 2 {ll.j— U 4 ) — a 4 («8— « 3 )— 5 ^°(«'7 — "2> I 

+ o 2 jo , M 13 — W8) + e G (»i2-"7)+g^(wii— «o) 

-ax«8-«3)-a 6 (« 7 -f2)- J a 8 ("6-«i)} 
- ^ a c ( w ,-« 2 )- 1 3 8 K-«o- J- 2 S 10 (m 5 - Wo )} 



=5 3 rc 3 « u + a 5 « 10 + -5'w s + J i «,„ + 5'« 9 + -a 9 w 8 



g^S+s 9 « 8 + ga'i« 7 j_fa3„ 6+S 5 M . + 1, 



+ ¥>u l0 + S 7 u a + -a 9 « 8 + a 7 ?/ 9 + a°« 8 + ^3'X 
o 

+ l($>u s + Z"u 7 + ] -S^u 6 \-h'u s + ^u i + ~Pu 3 
Vu 4 + cht 3 + ^Vht?) -l(Pu 3 + S»u 2 + U l3 «i 

+ J- 2 (a u «7+3 ,3 «6+ 3a i3 « 6 )-^^4+s 9 M 3+ ^a>s) 
- ^ (a <J «3 + S'X + 1 « ia «i) - ^ (a u «2 + s ,3 «i + 1 a i5 « )] 

= 53{a 3 ( M „-« C ) + 3a5( Wl0 - M5 )+ y « 7 («9-«4)+ y ^(«8-«3) 

+ J?S n ("7-%)+ |i« 13 («6-«l)+ g5«"(*-«o)} 

=5 4 {a 4 M 8 +4a 6 « 7 +c-8a 8 « 6 +6-4a ,0 « 5 +3'64a i2 M 4 +i-28a , ^3 

+ '272a lfl K 2 +-332a l8 tt, + -0016a 2 X} 
whence, dividing by 5 4 and transposing, 

(4) a 4 « 8 =A 4 -4a c « 7 -G'-8c 8 « 6 -G-4a l0 « i -3-64a 12 ! < 4 

— l-28a 14 « 3 --272a'«« <2 --032a l8 « 1 --0016a 2 o« ; 
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and this is in effect the fourth formula in Briggs's table, which we 
have inserted at art. 11, from page 29 of his Preface.* 

19. We shall not proceed any further with these calculations, 
which present no difficulty. Their law is evident, and they are 
eminently symmetrical, It must be admitted they are simple, 
although far from brief. All the operations are purely arithmetical, 
and depend on the proper use of the relations (B) and (C), applied 
to two direct consequences of the formulae (D) and (E). The 
minute resolutions required return with such symmetry that we 
can almost write down the results mechanically. 

We therefore think it highly probable that this was the course 
followed by Briggs. It does not presuppose the use of any know- 
ledge superior to that of a contemporary of Kepler and Harriot, 
whose works appeared about fifteen or twenty years before the more 
advanced era of Descartes and Fermat. 

But, as we have previously stated (arts. 9 and 12), the process 
appears to us rather curious than useful. There appears to be 
nothing to guide the computer as to the limit he should prescribe 
to himself in making his corrections ; and, moreover, the differ- 
ences of various orders having reference to different terms in the 
original series might give considerable trouble in practice. f 

20. On the contrary, nothing can be more lucid and convenient 
than the method published in 1670 by Mouton ; and although, in 
the last century, two different demonstrations of it have been given, 
we trust that the more simple and complete exposition we are about 
to give in the following section may not be considered superfluous. 

§ IV. 

21. " Given a series of quantities such that the differences of 
any given order are constant, to find any number whatever of 
intermediate values that shall follow the same law." 

Mouton was the first person who regarded interpolation from 
this point of view. The process published by him, after Regnauld, 
is a general one, but very long. He solves the problem by showing 
how to find by means of the differences of the original quantities, 
those of the interpolated terms ; these terms being then found from 
the differences by successive additions. Nothing can be simpler 
than this process. 

* M. Maurice's demonstration fails to exhibit the general law of this formula for 
£ 4 U8. This defect we hope to supply in our next number. — El>. J. I, A. 

+ It appears to us that M. Maurice does not do full justice to Briggs's method in 
consequence of his not observing that the differences o" 4 « 8 , ^ 6 "7, <Pt>6, S w uz .... are all 
"central differences" opposite the same value, «io. — Ed. J. I. A. 

VOL. XIV. C 
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The only difficulty then is to find these differences of the inter- 
polated terms, the number (h — 1) of the terms interpolated and 
the order (») of the constant differences in the given series, being 
any whatever. 

To solve it, Mouton and Regnauld had arrived at the following 
principle, probably suggested to them by the examination of several 
particular cases: — When the values of a function corresponding to 
equidistant values of the variable have their differences of any given 
order constant, then the values interpolated between them at equal 
intervals will also have their differences of the same order constant. 

This principle is true, but only in the case of rational and 
integral functions, as in Mouton's problem, and as we will now 
demonstrate. 

22. Let u be a function of x, such that 

u=ax m + (lx ! ' + yx'' + 

where the indices are all positive integers, and m>p>q> . . . 

Let us consider the relation u=ax m . Since m is the highest 
power of * in the value of u, the consideration of the other terms 
is useless for our purpose. Supposing Sx=h, we shall have 

Au=a{(x + h) m — x m }=amhx m - l + Ax m - 2 + Bx m - 3 + 

from which we deduce 

A?u=amhAx m - 1 + AAx m - 2 + 

and consequently 

A 2 u=amh(m—l)hx m - 2 + A'x m - i + 

In the same manner we shall find 

A 3 u=amh(m— \)h(m — 2)hx m - 3 +A"x m - ' + 

A"M=oBi(m-l)(m- 2). . .(m~n+l)h n x m - n +A i "- 1) x m -"- 1 + . . . 

Now, if m = n, since n>p>q> . . . , we see that if we had 
considered the complete value of U, all the terms of A"a after 
the first would necessarily have disappeared; so that we should 
have had simply 

A n u=n(n — l)(?t— 2) 2ah" 

for the value of the constant difference; for A n+, w = 0. 

It is also evident that if, taking §<# = 1, we represent by Su 
the difference (x+ l) m — x m , we should find in just the same way 

S"u=n(n— l)(n— 2) 2.1. a, 
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and we should consequently have 

6 U= —7— . 

h 

Thus the supposed principle is demonstrated, and the value of S"k 
is found. 

23. Assuming then that we have a series of quantities, w , «a, 
u ih , . . . , which give the differences Aw > A 2 w , A 3 m , . . . A n w > 
of which the last is constant, we see from what precedes that we 
should know also S n u , which relates to the series u 0) u u w 2 , • • • 
Uh, . . . and is, like A"k , constant. Having also m , if we could 
find hi , S 2 a , . . . B n ~ l u , we might by simple additions obtain 
the quantities m,, m 2 , w 3 , . . . , Ka-i> . . . 

To establish this, let us confine our attention, like Briggs and 
Mouton, to the case when n=A=5, and bear in mind that in 
general 

(»») «»+i=«m + Sw„. 

Thus, when u represents a number, it will be seen that each 
value is formed by the addition of the preceding one and of its 
difference; and as we know that u m — (1 +B) m u 0) it is clear that we 
shall have 

Consequently, Bu m will be expressed in terms of Su , and its suc- 
cessive differences; so that we shall have, m being any whole 
number, 

s ^ , 52 , m ( OT - 1 ) 5 , , m(m — l)(m— 2) w 
e«™= "m+i - w,„=K + m&i> + x 1 — - Pu -\ — i — yy Z* u o + ■•• 

It will therefore be sufficient, as already stated, to know how to 
express the differences Su , <$ 2 u , § 3 w , S 4 m , in terms of the differ- 
ences A (S 5 m being a constant and its value known), in order to 
be able to calculate all the interpolated numbers. 

Assume for a moment that we know these four quantities, Su , 
. . . B*u , and we shall see how easily we may obtain the inter- 
polated numbers, by virtue of an obvious consequence of the relation 
(m) which gives 

(«) a l « m+1 =s'« m -i-a'+' Mm . 

24. Let us consider the following table, the law of which is 
evident : — 



c 2 
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Ora 


Interpolation. [Oct. 


«0i 










M, 


S«o> 








U 2 


Su { 


a 2 «< » 






"3 


Su 2 


2 2 W, 


2 3 «/ , 




«4 


Su 3 


S 2 « 2 


5 3 i/! 


ax, 


"5 


Su 4 


S% 3 


Pu 2 


£ 4 Wi 3 5 w (=constant, by hypothesis) 


M 6 


S«5 


Pu 4 


2 3 « 3 


2 4 M 2 <5 5 «0, 


M 7 


2«6 


a 2 " 5 


a 3 u 4 


2 4 « 3 S 5 «0, 



Every term in this table will be known, if we know only those 
at the end of each horizontal line. For, by virtue of the relation 
(n), in accordance with which it has been formed, each term being 
equal to the sum of the two which stand respectively above it and 
on the right of it, we see that 

gw, = gw +S 2 u 

u 2 = u x +&l\ 

<S 2 w t = Puq+Puq 
■ $u 2 = S«i +^ 2 «i 

U 3 = « 2 + &h 

Put = S 3 u + B*u 
lht 2 = 8 2 «i + t> 3 «i 
lu 3 = Su 2 +Pu 2 

«'4 = u 3 +2»3 

and so on ad infinitum. It is furthermore easy to see how we 
must proceed if the constant differences are of either a higher or a 
lower order than &u . 

Thus, each of the interpolated values will be formed by adding 
to the previous one its difference, which itself will be obtained by 
simple, regular and successive additions, and will involve §u , § 2 w , 
B 3 u ... up to the order for which these differences are constant. 

25. We have lastly to find Su , S 2 m , . . . §"«„, in terms of Am , 
A 2 w , . . . , A"w , supposing that A"m is constant. 

Let there be two series of values of the function, in one of 
which every value is known, but in the other only every hth value 
is known, namely : — 

w , u x , a 2 , . . . . u h , u h+ i, . . . u ih , u. 2h+1 , .... u 3h u, h , .... 

It is clear that if v =u , t'! = M A , we shall have generally, 
v x =n xh . 

Denoting by § the differences of u, and by A the differences of 
v, we have 
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«i=«o+<$«<» i>i=t> + Ar =Wo+A»>o, 
whence Av =v 1 — u =u h — ie . 

We therefore have, since m a =(1 + 8)*m , 

Ai> ={(l+a)*-l}«o- 
Let us confine ourselves, for simplicity, to the case of A =5. 
We shall see however that the course of procedure is quite general. 
Let us also put (1 -+- S) 5 = /3. We shall then easily get 

v l =u„=fiu ' A« 1 -Ar =A 2 fo=(/3— l)%o=A% 

ai A»,=0(/3— 1)k„ .„ ,,, 1V > 

»2=«2*=/3 2 «o . n%n ,s A 2 t',=(/3-l) 8 « 

, 5=MM =^ A- 4 -^-!H 



A^ 2 = W -l)3 ao AV = ^/3-l)S 



We have thus shown that 

Ar =Aw =(/3-l)« ={l+a) 5 -l}M 
A 2 f =A 2 a =(/3-l) 2 Wo ={l+5)5-l}2 M() 
A»r =A»«o=(/3-l)»ii ={l+3)»-l}«K 
AS=AX=(/3-l) 4 %={l+S) 5 -l} 4 «o 
A 5 f =A s M =(^-l)X={l + e) 5 -l} 5 «o 

and it is clear that, whatever integers m and h may be, we shall 

have 

A™M ={(l + a) 4 -l}'"«o; 

whence expanding { (1 + S)*— 11 , we have 

(P) [, S + ^W^<^3 + ....p 

=/rs ra + As m+i +Ba'"+ 2 +ca m+3 + 

If then we can by any means determine the numerical values 
of the coefficients A, B, C, . . . , we shall only have to multiply the 
second member of (P) by m , and we shall have A m u in terms of 
S m u , S m+I u > • • • As to the determination of those coefficients, we 
shall see that the law of their formation is very simple ; and we 
confine ourselves here to remarking that if the differences S m be 
constant, the differences of higher orders will vanish ; so that we 
again get A m M =A m g m M . 
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26. But the problem is not to find the differences A in terms of 
the differences &; it is the converse problem we have to solve. 

For that purpose, since v t = (l + A)» = m&=(1 + §) h u , we get 

(l+c)« =(l+A)%; 
or, since u = v , 

hi ={{l+Ay-l}u . 
Reasoning as in art. 25, we shall have 

3'%={(1+A)*-1}>, 
Hence, when m=l, we have 

r , , i ,, fA 1(1— A) A 2 1(1— h)(\-2h) A 3 ) 

But, if we take A = 5, suppose A 5 « constant, and for brevity put k 
for t" , we can very easily calculate the numerical values of the 
coefficients of k 2 . . . k 5 , and we thus get 

(1 ) hi =(k+aP-\- blfi + ck 4 + <W)i< . 

But S 2 « will depend on the square of the value of B thus found ; 
and since we need not consider higher powers of k than A: 5 , putting 
k + ak 2 = a, bk 3 + ck 4 + dk 5 —(i, we see at once that in the 
expansion of £ 2 w = (a + j3) 2 w , we may confine our attention to the 
following 

g%=(a 2 + 2c43K, 
which gives 

(2) 2 2 « = {k 2 + 2aP+(a 2 +2b)k 4 + 2(ab+c)k*}ii . 

In the same way 

Pu =(a+l3fu 
reduces to 

c 3 e< =(a 3 +3a 2 /3)%, 

whence we obtain 

(3) Pu = {£> + 3«A- 4 + 3(« 2 + b)P}u . 
In the same way, B 4 u — a 4 u leads to 

(4) a% =(^ + 4«A> 0) 
and lastly 

(5) a^ =p« =^-\ 

as we have already more than once seen. 
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27. The actual calculation of all these coefficients is as simple 
as it is rapid. Thus, when A=5, we have 



whence 



jiz^ = _21= e , ci^ = +79-8=rf; 



2a=— 4, a 2 + 2A= + 16, 2(rt5 + c)=-66, 
3a= — 6, 3(« 2 + J)= + 30, 4a=-8. 



Thus, putting for A its value v , we shall have 

o 

(i) a Mo= ^_ 2 ^.o +6 ^_o_ 21 ^.o +79 . 8 ^o 
(2) ax= ^o_ 4 ^o +16 ^_ 66 ^o 

0) *".- ^- 6^ + 30^-0 



(4) 3<« = 

(5) 2 5 « = 



A<</ A% 

AX 

5 5 



Consequently, since we are supposed to know the numerical values 
of Aw » • • • A 5 M ,we shall very easily find the numerical values of 
Su , S 2 u 0> . . . 8 5 « ; and. all that is necessary to form the table in 
art. 24; so that, the origin u being given, all the interpolated 
terms will be found by very simple successive additions. 

We further notice that while it will never be necessary to con- 
sider the differences of a higher order than the seventh or eighth, 
even in these extreme cases, by following the same process, we shall 
readily be able to obtain all the numbers necessary for the process 
of interpolation. As to the value of h, however large it may be, 
no difficulty can arise. 

28. It will be interesting to compare the five formulae above 
given with those of Briggs, demonstrated in arts. 15 to 18. 

But since Briggs has considered S'"u n we must remember that 
since 

w„=(l + g)"« , 
we have also 

(a) tfx=S" , (i+c)"M . 

We must also bear in mind that having proved (art. 25) that 
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the lowest power of k or =- in the value of S'"m is m, it will be 

o 

useless to carry the expansion of the formula (§) beyond § s w ; for 

&u 0) B 7 u , &c, will involve no lower order of differences than A 6 w„, 

A 7 w , &c, respectively; and these are all zero, because A 5 m is 

constant. 

Then, taking the formula (8), the value (1) of art. 15, 



<VW 2 = A 1 —c 3 ti\ — -= £ 5 «o 
5 



will become 



o 
whence 

(1') Bu = ^2 _ 2J% -2Pu -S*u - J * 5 » . . 

In the same way, the value (2) of art. 16, as far as it need be 
taken, will be 

S% 4 =A 2 — 2a 4 « 3 

or Pn +4dhi + 6S"u + 4S 5 </ = ^-° -2« 4 i« — CS-''k 0> 

whence 

(2') 8*« = ^2— 4a»« -8S«if — ioa»Mo. 

The value (3) of art. 17 will similarly give 

t 3 !f 6 =A 3 — 3<S 5 ?< 6 

or S 3 « + 6S 4 « + 153 5 w = ^-° - 32 V , 

whence 

(3') Pu = ^ -62^< -182 5 « . 

The value (4) of art. 18 will give in the same way 

3 4 «o=A 4 , 

or a4„ o + 8 a5 Mo= ^_"'> ) 

o 4 
whence 

(4') Vu =^-Wu . 

Lastly, Briggs having shown (art. 11) that the last difference 
needs no correction, it is clear that he would take 

AS> 
5 5 ' 



(5") 3*ir =A»= r5 



1867.] On Interpolation. 25 

We shall find then, finally, by successive substitution 



(4") 


S4„ A4 "° 

«S>- 54 • 


„AS> 
8 5» 


(3") 


S3 AS) 


"6^+30^ 


(2") 


*s,=- 5 - r - 


.4^+16^-66^ 
5 3 5 4 5 5 


(1") 


= ~z~ ' 


-^ +6^-21^+79-8^ 



and we see that tbese values (1") . . . (5") coincide exactly with 
the values (1) . . . (5) of art. 27; so that the methods of Briggs 
and Mouton, notwithstanding their great dissimilarity, lead to 
identically the same results. 

§v. 

29. Hitherto we have supposed the given values to be equi- 
distant ; but this is not always the case. 

Thus, for example, in observing certain phenomena it some- 
times happens that the intervals between the observations are not 
equal. In that case, the time being taken for the variable, and 
the results of observation giving the values of the function u x , those 
values will correspond to values of the variable, which do not pro- 
ceed by a constant difference ; and the problem consists in finding 
a general expression for u x in terms of the given values and one 
value of x, corresponding to any time we choose. 

Let the known values of u x , corresponding to the values of x, 
0, p, q, r, s . . . , be u , u p , u q , u r , u s . . . ; then we have to find a 
general expression for u x , whatever be the value of x. 

Reasoning as before, we may assume 

U;=u + Bx + Cx(x— p)-\-T)x{x — p)(x— q) + Ex(x—p)(x—q)(x—r) + ... 

and then denoting by A, A', A", A'" . . . the differences u p —u 0) 
u q —u , u r —u , ... we have 

«,— «o=A=Bj», u,—u =A'=Bq + Gq(q—p), 
u r —ti =A"=Br + Gr(/—p) + Dr(r—p)(r—q). . . . 

Whence we shall find 
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(«) 



A 

P 

A' 



?(?—i>) P{P~'l) 

A" A' 
D = "7 U 7 + -7 r, ; + 



Kr—PXr—S) l{<l—P)(<i — r ) P{p — t]){p—r) 

A"' A" A' 

+ ~, T7 CT -x + ' 



s(s—p)(s-q)(s—r) r(r—p)(,— q){i—s) q(q—p)(q->)(q—s) 

A 



+ 



P(P—<i)(P—>)(P—«) 
&c.= &c. 

the law of these successive values being simple and obvious. This 
symmetry, as well as the recurrence of the positive sign which 
unites all the terms, results from the fact that in the several factors 
of each denominator we have throughout written that letter first, 
which is itself one of the factors. Thus, for example, wc have 
written p(q — p) under the form — p(p — q) '■ an d similarly with the 
others. 

Having thus obtained the values of the coefficients B, C, D . . . 
we shall have for the required expression, collecting the terms 
involving A, A', A" . . . respectively, 

C x x(x—p) x(x — p)(x — q) x(x — p)(x — q)(x — r) ) 

v '~ "° \p p{p—q) p{p—q)(p—r) p{p—g)(p—r)(p—s) ""} 

x(x—p) x(x—p)(x — q) x(x—p)(x — q)(x—r) ) 

?(?— p) q(q—p)Kq— r ) q(q—pX9— r )(q— s ) ' "1 

| x(x— p) (x-q) x(x—p)(x—q)(x—r) } ^„ 

\r{r— p){r— q) r(r—p){r-q)(r—s) '"J 
( x(x— p)(x— q)(x— r) ) ^ 

1 s(s—p)(s—q)(s—r) '"} 

30. Let now X = x(x — p) {x — q ) (x — r)(x — s) . . . , and denote 

by X w . Then we should have similarly 

x — w J J 

V P=p(p—q)(p- r )(p— s ) 

Q,=9(i~p)(i— r )(<i- s ) — 

R r =»'(r — p)(r — q)( r — s ) ■ ■ • • 
S]=s(s—p)(.s—q)(s — r). . . . 
&c. = &c. 

If we next consider the coefficients of A, A', ... , in the above 
formula (a) for u x , we shall notice that the several terms in them 
can be grouped according to a regular law. 
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Thus, the first being, for example, 

x(x — p). . . .{x — f) 



v(v — p) .... (v — f) 
the sum of the first and second will be 



orV, 



X V\ ~X V 

1 H ! — V - =Vj suppose. 

v—y, v—y 

The sum of the first, second, and third will be 

V, ( 1 -| ) = V, - — = V 2 suppose, 

V v — zj v — 2 

and so on, where we see the simple factor x—v has disappeared 
entirely. Treating each of the coefficients of A, A', A", ... in this 
way, we can easily see that the value of u x can be accurately written 
in the concise form 

(«') u,=« + £ A + £ A'+ £ A" + f A'"+ .... 

i p % iv r b, 

In this formula, each of the coefficients of A, A', A", . . . (which 
we will denote by a, /3, y, 8, . . . for brevity) is easily calculated by 
logarithms ; and the number of the coefficients will be the same as 
that of the values p, q, r, . . . 

31. This expression for u x has been readily obtained, because 
we have considered only the differences u p —u , u 9 —u , u r —u , 
&c. ; but since these differences of necessity increase, the expres- 
sion is not convergent, and we shall find a more convenient form 
by introducing instead of A, A', A," ... the differences Aw , AWp, 
Am ? , . . . between the successive values of u x , which will not 
necessarily increase. 

For this purpose, we notice that 

A —u p — w =A«<o 

A' = u ? — w =m,— u p -\-u p — u —Au p + Au 

&" = Ur—Vo= u r—u <l +uj—u p +ii p —u =Au 7 + Au p + Au 

and so on. 

If, then, we make these substitutions in the above value of u x , 
denoting as before the coefficients by a, /3, y, . . . , the formula (a') 
will take the preferable form 

(a") u,=u + (a + (} + y + S+ )Au + ((3 + y+S+ )Au p 

+ (y + S+ . . . .)Au,+($+ . . ..)A« r + .... 

32. But we can very easily obtain a symmetrical and more 
convergent formula. In fact, by the last article 
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(1) 



whence we obtain 



A = Aw 

A' = Au p + Au 

A" = A?/, + A' 

A'" = Au r + A" 
A IV = A?>, + A'" 
&'C. = &C. 



(2) 



A« = A, 
Au p = A'— A, 

Au, = A"— A', 
Au, = A'"— A", 
Au s = A' v — A'", 



Au p — Au = A' — 2A = A 2 ?/ , 
An,— Au p = A" —2 A' +A = Ahi p , 
Au r -An, = A'"— 2A" +A' = A 2 u,, 
Au s —Au r = A ,T — 2A'" + A"= A 2 u r , 



(3) 



A 2 u p — A 2 u = A" — 3A' +3A =A 3 u , 
A 2 u,—A 2 u p = A'"— 3A" +3A— A=A 3 « P , 
A 2 u r -Aht, = A ,V -3A"' + 3A"-A'=AX, 



(• A 3 « ? -A 3 w =A'"-4A" + 6A' -4A=A 4 « , 
(4) { A 3 ?( 5 — A 3 ^=A ,V — 4A'" + 6A"-4A=A 4 m,,, 



1 

(5) { 



A\-A i u =A"— 5A'" + 10A"-10A'+5A=AV , 



From these we readily conclude that 

A = Au , 

A' = A 2 « +2Aw , 

A" = A 3 M + 3A 2 « +3Aw , 

A'" = A 4 « + 4A 3 « + 6A 2 « + 4Aw , 

A" = A 5 m + 5A 4 « +10A 3 « + 10A 2 « +5At/ , 

So that, substituting these values in (a'), it will become 

(a'") % = Wo + [ a +2/3 + 3y + 4S + 5£+ . . . .].A« 
+[/3 + 3y + 68+10 £ + . . . .].A 2 « 
+ [ y + 43+10£ + ...].A 3 « + [a + 5£ + ...].A 4 K 
+ [ £+ .... ].A S « +..-., 

a perfectly regular series, of which the law is evident, the differences 
of u in general converging rapidly, while the computation demands 
only, so to speak, that of the numbers a, /3, y, 8, e, . . . which is 
very easy. 
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As to the values of A« , A 2 w , . . . A 5 « , . . . the equations 
(1). (2), ... (5) give each of them in terms of the immediate 
differences A, A', . . . A ly . . . 

33. This was not the form in which the first solution of our 
problem was given by Newton in Lemma V. of the third book of 
his Principia. Laplace, in his Mecanique Celeste, has given a 
solution virtually the same as Newton's, but in a rather different 
form, which we will now demonstrate. 

We will, in the first instance, establish an Algebraical Lemma, 
which for the sake of symmetry we shall find it convenient to use 
in certain subtractions, and which may prove useful on other 
occasions. 

" Given any number, m, of binomial factors, as x — a or a, x — b 
or j3, x — c or y, &c. . . . ; and the same number involving y 
instead of x, but with the same second terms, viz., y— a or A, y — b 
or B, y — c or C . . . , the difference oj3y . . . —ABC . . . will be 
divisible by x—y; and the quotient will have two equivalent sym- 
metrical forms." 

1°. We see at once that the difference will be a polynomial of 
the form 

^-jr-p(^- I -/- I )+Q(^- i -ir-*)-R(^- , -jr-'). • ±w.(*-y), 

and that this polynomial is divisible by x—y ; 

2°. If we suppose x—y = z, it is easy to see that the differences 
a— A, /3— B, y — C, .... are each equal to z, and that we shall 
thus have, ad libitum, 

either (1) a=A + z, /3=B + z, y=C + z , 

or (2) A=a — z, B=/3— z, C=y— z , 

Thus, according as we employ the relations (1) or (2), and take 
m equal to 2, 3, 4, ... , we shall have the following as the differ- 
ences of the products under consideration : — 

Successive According to the According to the 

values of relations (1). relations (J). 

a/3 -AB = (B + a> = wz = (/3 +A>, 

a/3y —ABC = (Cir + a/3)z = 7r,z = ( y! r +AB)s, 
a/3 y a -ABCD =(D*-, + a/3y)z = 7r a z = (tV, + ABC)z, 
a/3yg£— ABCDE=(Ejt 2 + a(iyS)z= 7r 3 z = («r 2 + ABCD)z, 

We have thus obtained two sets of values for the quotients ir, 
iT], 7r 2 , tt 3 , . . . the expansion and comparison of which will give 
the following symmetrical relations : — 
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B +a =7T =/3 +A, 
BC+ Ca +aft =7r,=/3y + yA +AB, 
BCD +CDa +Da/3 +ufiy =^ 2 =/3yS + ySA + SAB + ABC, 
BCDE + CDEa + DEa/3 + Eafiy + a/3yS=7r 3 =/3 y a £ + y£ £ A + g £ AB 

+ £ ABC + ABCD, 

and it may be observed that either of these two series may be 
derived from the other, if in that which vvc assume as given, we 
write a, j3, y, . . . for A, B, C, . . . , respectively, and A, B, C, . . . 
for a, /3, y, . . . , respectively. 

31. So much being premised, let n, p, q, r, s, . . . . be values 
of x to which correspond the following given values of u x , viz., 
u, u, u", u"', u". ... It is clear that we may assume 

u x =u + B(x — n) + C(x — n){x — p) + T>{x — n){x — p){x — q) 
+ E(x—n)(x—p){x—q)(x—r)+ , 

where the coefficients B, C, D, E, . . . are to be determined by 
means of the given values of u x . 

If we make n — in the above value of u x , it becomes identical 
with the one in art. 29; but we have preferred the abore form, 
because it leads to more symmetrical results. 

Now putting x equal to n, p, q . . . successively, we have 
u =u, 

u' =« + B(/) — »), 
u" =zu + B(q — n) + C(q — n){q— p), 

u'"=u-\- B(» — w) + C(> — ii)[r — p) + D() — n)(r—p)(r — q\ 
u lw =u + B(s — n) + C(s—n)(s—p) + T>(s — ?i)(s—p)(s—q) 

+ E(s—n)(s—p)(s—q)(s—r), 

> 

From these equations we can deduce by means of the lemma 
in art. 33, 

u — u =B(jB — n), 

u" —u = B( q— ])) + C(q—n)(q—p), 

«"'-»" =B(r— q) + C(r—q).[(q—ti) + (r—p)] + D( r —n)(r—p)(r-q), 

u "—u'" = B(s — r) + C(s—r).\_(r—p) + (*—«)] 

+ D(*— r)[(r— n)(r—p) + (r-n)(s-q) + (s-p)(s~q)2 

+ E(s-n){s-p)(s-q)(s-r). 

Now dividing and using the symbol S for brevity, we obtain 

91 —. tl 

-.Su =B, 



p — n 
u' — u 
q—p 



=hi =B + C(q—n), 
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tu i' 
« — u 



r-q 

It" — II 



=Eu"=B + C(q-n) + C(r-p) + D{ r ~n)(r-p), 



-.Eu'"=Ji-\-C(r—p) + C(s — >i)+'D(r—n)(r—p) + D(r-n)(s-q) 
+ T)(s-p)(s -q) + E(s - n)(s -p)(s - q). 



A similar operation will give also 

=E 2 u =C, 



on — cu 



q-i 
En" — En' 

i — p 
En"' — Eu" 
s-q 



=E 2 u =C+D(V-»), 

=E 2 u"=C + !)(>■ -n) + D(s— p) + E(s— «)(*— p), 



and in the same way we shall find 
Ehif—E^/i 



=E 3 u —D, 



r — n 

c' 2 «" — E 2 u 
s — p 



:oV=D + E(s- n), 



3 u' — c 3 u 

=3 4 m=E. 



The other coefficients, if required, would be found in exactly the 
same way. 

Substituting now the values of the four first coefficients, as 
just found, we have for the value of u x 

u x =u-\-(x — n).Su + (x — n)(x — p).c 2 u + (x — n){x — p)(x — q).E 3 u 
+ (x — n)(x — p)(x — q)(x — r).E 4 u-\-&c; 

The process of computation by means of this formula is very 
easy, the numerical values of the quantities Sw, S 2 «, . . . being 
known, and being in most instances sufficiently small to allow of 
our confining ourselves to a few of the first terms of the series ; a 
point easily decided in any case. This is the form adopted by 
Laplace, Mecanique Celeste, torn, i., p. 200. 

35. Notwithstanding the dissimilarity between this solution 
and ours, it may be shown that they agree exactly; and we will 
conclude this paper by proving that such is the case. 

If we substitute for Zu, S 2 w, S 3 w, . . . their values deduced from 
the mode of their formation, as explained above, we shall find 
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u' u 

Su = 1- 



p — n n — p 

Pu= + 



Shi= 



(q-n)(q-p) (p-n){p-q) (n—p)(n—q)' 

u" u" t u' 

<,r—n){r—p){r—q) (q—n)(q—p)(q—r) (p-n)(p — q)(p—r) 

u 



+ 



(n-p)(n—q)(n—r)' 



Nothing can be more symmetrical than these expressions. 

Now in order to compare this solution with ours, we must, as 
already observed, make n = 0; and we shall then have, for example, 
u'" u" u u 

* jW ~ ~r(r-p)(r—q) q{q—p){q—r) p(p-q)(p — r)~ pqr' 

which, for brevity, we will write 

u" u" u' u 
= ~R + Q + ¥~p~q~^' 
But in art. 29 we found the following value for D, which ought 
to be equal to the preceding, since we have put n=0, 
„ A" A' A 
D= K+Q+P' 
„ u" — u u" — u u — u 

We must therefore have 

1 1, 1_J_ 
K + Q + P~pqr' 

which is seen to be the case, when for R, Q, P, we substitute their 
values. 



APPENDIX. 

On a casn that may occur in astronomical interpolations. 

I. Wlicn a series of astronomical observations is made, the element 
observed may be considered as a function of the time, since it varies there- 
with. Let t denote the time. We may then, as in the calculus of finite 
differences, represent the value of the element undei consideration, whether 
observed or computed, by u, ; the time t being considered as the independent 
variable. 

Suppose that we know n values of the element, or of the function u lt 
viz., a, a', a", . . . corresponding to the times 0, &, 6", . . . which are 
supposed to be successive and equidistant, and that we wish to deduce a 
general expression for u t that shall be true for any value of t not too remote 
from the period under consideration. The required solution will be 
possible, and indeed very simple, if the successive differences of the given 
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values become either constant, or small enough to admit of our neglecting 
all beyond a certain order. 

Suppose, for example, that this is the case with the fourth differences, 
so that we may neglect the differences of the fifth and higher orders: it 
will be sufficient to consider five values of n„ or to have »=5. 

Also let 

ff-8=6"-0'=ff"— ff'=0"-ff"=a. 

The common difference of the increasing arithmetical progression 
6, &, . . . 6", will then be a, and it is clear that we may assume 

u t =A + 'B(t-6) + C(t-8)(t-d') + D(t-6)((-#)(t-ff') 

+ E(t-6)(t-6'){t-6")(t-0"'). 

The given values of u, will then lead to the following equations: — 
a = A, 
a! = o-(-B.a, 
a" = a+B.2a + 0.2a 2 , 
a'" = «+B.4a+C.G« 2 + D.6a 3 , 
a" = a+B.4a+C.12a 2 + D.24u 3 +E.24a\ 

and we may thence conclude almost directly, from the known values of the 
differences of any order in terms of given numbers, that 

A=a; B=-.Sa; 0=-^-^; D= — \-. -Pa; K = T ^-J-j— 4 2 4 «. 
a 1.2. a 2 1.2. 3. a 3 1.2. 3. 4. a 4 

We shall thus have for the solution sought — 

m ^'- e * , (t- eyt-tf) {t -0)(t-ff)(f- 6") 
( A) tt _ fl+ _. 3a+ ___ 5 — a 2 «+ -— a»« 

(< _fl )(< _ 8 ') (< _r)(<- r) 
+ IZSZ^ ca - 

It is also easy from the symmetry to perceive the form this formula 
would take if n were larger, and we were called upon to consider differences 
of a higher order than the fourth ; and we see readily what terms should 
then be added to the preceding. 

II. Now if we desire to reckon the time t from the middle of the period 
under consideration, so that 0"=O, and to take the time between two con- 
secutive observations as the unit of time, we must, in the formula (h), put 

«=1, 0=-2, ff=-l, fl"=0, 0"'= + l. 

That formula will then give 

(Ji) t, l=a .+ —.Za-\ — -b 2 a+ YJ^ ,Pa 

, (t+2)((+l)f( t-l) 
1.2.3.4 
where, as before, we have only to do with a and its differences, as found 
from the five given values of u,. 

III. Stirling (followed by Lacroix, tome iii., pp. 27-8) has given the 
formula for a similar case under a different form, which he has only been 
able to generalize by induction. M. Bessel, who has employed the formula 

VOL. XIV. D 
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given by Stirling, appears to have demonstrated it in the following manner 
for the case where «=5. 

Let u^A+M+Cfi+DP + El 4 . 

Then making I successively equal to — 2, — 1, 0, -|-1 and +2, we shall 
have five equations which we may denote respectively by ( — 2), ( — 1), 
(0), (+1) and (+2). 

Then (+1) + (_1) gives ?-±- -«" = C + E; 

(+2) + (-2) „ £±f _«" = 4C+1GK; 

(0) „ a" = A; 

(+l)-(-l) „ 'LTJL =B+I); 

( + 2)-(-2) „ "-f- =2B+8D. 

We hence obtain, by means of the well known values of the various 
orders of differences in terms of the given numbers a, a', . . . a", the fol- 
lowing values : — 

" ,» hi -\-lu" Pa-\-Pd l 2 d e 4 a 

A=«, 1J=— 2 - -— , 0=----, 

2.6 ' 24 

Substituting these values, and rearranging the terms with reference to the 
powers of t, we get 

r,V + cV'-] / t* rPa + Fa"l /(7 2 -l) 

But this formula is subject to the inconvenience of involving differences 
of several of the given values; while the formula (A') contains only the 
several differences of the first value a. It also appears much less 
symmetrical, and less easy to remember : and, above all, it is more difficult 
of rigorous generalization than the formula! (//) and (//). 

Notwithstanding the apparent difference of tlie formula? (//) and (y), it 
is easy to prove that they only give the same result under different forms. 

Another case of astronomical interpolation. 
IV. In astronomy the third (if not the second) differences are generally 
constant, or, at all events, small enough to permit of our neglecting the 
fourth differences; and the formula (A) then reduces to 

a 1.2.(1'' 1.2.3.a J 

in which d, 0', 6", ()'", being equidistant values of /, corresponding to the 
known values a, a', a", a'", of n„ we have 

a = 0' - 6=0" - &=&'' - 6". 
Now, if the third differences are constant, we take from the tables, or 
we know from the observations, only the four values a, a, a", a", of which 
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the former two will correspond to times that precede, and the others to 
times that follow the epoch t for which we are computing, and which must 
consequently fall between ff and 0". It will therefore now be convenient 
to make 0'=0, whence 

0= — a, and 6"= + a. 
The formula then gives 

t+* , U+a)t „„ t(l 2 —a 2 ) „ 
a 1.2. a'' 1.2.0. a" 

or, 

a l.z.a' 4 l.z.oa* 

/ < ,,.<(<-") l,-.. , «.„ -, t(t~a)[~t+a 11 „ 

L. a a' 4 1.2. 3.0^ 

This formula (/) coincides exactly with formula (2) of p. 99 of 
L' Astro nomie pratique of M. Francceur, when we make 

a=12\ whence -a=6\ 

and put A'=Ja', and ^>= - (3V + S 2 a). 

It may also be made to coincide with the formula (3) at the top of 
p. 100 of the same work. In fact, the correction x or it, — a, to be applied 
to the value a, becomes 

(») «=i[y_*+ 1.3» a ] + £l[ _l.g3«J + '±.l.Pa, 

which is of the following form — 

I fl t 3 

*=A — +B-— +C-, . 
a a-* a* 

If now we pursue the reverse process and start with this formula, we 

shall find the values of A, B, and C, by making t equal to — a, a, 2a 

successively in the equation 

t t" 1 t 3 

m,=o' + A--+B.— +C- — . 
a a' a 3 

We thus get 

a = a'-A + B — C, 
a"=a' + A+B + C, 
a"=d + 2A + 4B + 8C. 
The two first of these equations give 

«"+a-2ra'=2B, 

whence B= - • Pa= - (Pa + S 2 o) Pa. 

2 4 4 

Also a"— a=2(A + C). 

But, from the third, a'"— 4B— o'=2A + 8C; 

d 2 
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whence by subtraction, 

a'"— a— (a— a)— 4B=6C, 

and ha" — la — 2.3 2 a=6C; 

or, since la" — la = I 1 a -\- l 2 a = 2 . l 2 a + l 3 a, 

we have C= . 

l.l.O 

Finally, as the second equation gives 

A=a" — d — B — C, 

we shall have A=Sa'— <f>+ — • S 3 «, 

and substituting these values of A, B, C, we get the formula (w). 

It may be noticed that in using the formula (/), if we neglect altogether 

the last term which involves Pa, the second differences would still be partly 

corrected by the differences of the third order. In fact, the formula would 

then give 

, t „ , t(t— a) 

u=a ^ la H — '•$. 

a a 1 

But, 

so that we shall have 

It is from this last term that M. Mathieu has computed a table inserted 
in the Coniiaissancc des Temps. Supposing a=12 4 , and that the motion 
of the star is uniform during those 12 hours, the two first terms 

d + —j . da' are given; and we then apply the correction resulting from 

the last term by means of the table in question, which is computed for 
values of t differing by 10 minutes. 

Generally speaking, in all these formula!, we give to a such integral 
values as we mav require, making it most frequently equal to 12*, 6*, 3*, 
or 1*. 
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XlAVING thus found the required formula for the value of 

a^», viz.:— 



Tr-^ 2m 12m' 

1 D", 



r(l— mr)(l— 2wr) 



